Light scattering and localization in an ultracold and dense atomic system 
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The quantum optical response ol high density ultracold atomic systems is critical to a wide range of 
fundamentally and technically important physical processes. These include quantum image storage, 
optically based quantum repeaters and ultracold molecule formation. We present here a microscopic 
analysis of the light scattering on such a system, and we compare it with a corresponding description 
based on macroscopic Maxwell theory. Results are discussed in the context of the spectral resonance 
structure, time-dependent response, and the light localization problem. 
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Early studies of formation and dynamics of ultracold 
atomic samples were largely focused on obtaining the 
high density and low temperatures necessary to attain 
Bose-Einstein Condensation in atomic gas samples 
At the same time, a large number of other research ar- 
eas emerged from studies of ultracold atomic gas sam- 
ples. Among these, there are several where physical 
processes are importantly modified at high atomic den- 
sity. These include efforts in ultracold molecule forma- 
tion image storage in high-optical depth samples [1], 
and hght storage and manipulation for possible atomic- 
physics based quantum repeaters 4] and other quantum 
information appHcations. Another little explored area 
which may have significant impact in a range of scien- 
tific or technical areas is study of quantum optical pro- 
cesses at high density no, where n^X^ ^ 1 [5]. Under 
these conditions, disorder-dependent phenomena such as 
an Anderson- type [6] localization^andom lasing [3] , and 
other mesoscopic phenomena I^, [13] may appear. In 
these cases, the quantum optical response of the medium 
is of critical importance to dressing and probing the phys- 
ical processes involved. 

In this letter we are concerned with the collective dy- 
namics of light and a sample of ultracold and high-density 
atoms. Under the conditions that n^X^ > 1 the inter- 
action of atomic dipoles via longitudinal and transverse 
electric fields strongly interferes with the scattering pro- 
cess. The atomic dipoles cannot be considered as inde- 
pendent secondary sources of the scattered waves freely 
propagating through the sample, as usually takes place 
for an optically dilute system. Here we present theo- 
retical discussion of the scattering process based on two 
different and complementary approaches. In the first, 
we use a self-consistent description of the atomic sam- 
ple in the spirit of the Debye-Mie model for a macro- 
scopic spherical scatterer consisting of a homogeneous 
and dense configuration of atomic dipoles. We apply 
a self-consistent calculation of the macroscopic permit- 
tivity based on its relevant statistical expression given 



by the Kubo formula ll|. In this model the scattered 
light is restored as secondary waves created by randomly 
distributed atomic dipoles, whose reradiation should be 
properly expanded in the spherical modes of the Debye- 
Mie problem. In the second approach we make exact 
numerical analysis of the completely quantum-posed de- 
scription of the single photon scattering problem. The 
atoms, the incoming photon and the field modes partici- 
pating in the process are considered in a second quantized 
formalism. This analysis gives us the exact value of the 
scattering amplitude and in particular shows the correct 
description of the cross section in those spectral domains 
where the Debye-Mie model fails or is not irreproachably 
applicable. We show the presence of collective micro- 
cavity structure, created by the atomic system, where 
a number of resonant modes have a sub-radiant nature, 
and discuss their physical substance in the context of 
light scattering and diffusion. The results have impor- 
tant implications for studies of light localization, and for 
manipulation of single and few photon states in ultracold 
atomic gases. 

Light transport through a disordered system of atomic 
dipoles can be discussed in the Maxwell theory as a 
combination of the process of coherent propagation, de- 
scribed by configuration averaged Maxwell equations, 
and incoherent losses initiated by scattering from meso- 
scopically scaled permittivity fluctuations. In such a for- 
malism one can visualize the light state as a local plane 
wave existing on a spatial scale ^ X and coherently in- 
teracting with the system of atomic dipoles located in 
a relevant mesoscopic volume. The dipoles interacting 
with the field are locally indistinguishable and the in- 
teraction process is assumed to be well approximated 
by the cooperative coherent response [lH for the sam- 
ple susceptibility x{'^) frequency uj. Then the per- 
mittivity = 1 -l- 47rx(tt') can be extracted via a self- 
consistent approximation of the collective dipole dynam- 
ics driven by the transverse electric field and modified 
by the Lorentz-Lorenz effect coming from interference of 
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proximal atomic dipoles. If incoherent losses have only a 
radiative nature and atoms equally populate the ground 
state Zeeman sublevels we obtain a self-consistent expres- 
sion for e{Lu) for an infinite sample 
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This equation can be analytically solved and e{uj) ex- 
pressed in terms of well defined external parameters, such 
as the natural decay rate 7 and tiq. The reduced matrix 
element dpoF for a given dipole transition can be also ex- 
pressed by 7. We consider here a closed transition, such 
that for hyperfine levels only the selected lower and up- 
per states with angular momenta Fq and F respectively 
can be radiatively. The above result can be applied to 
calculation of the scattering cross section of light by a 
spherical sample of large radius a (a ^ A) in the stan- 
dard formalism of the Debye-Mie problem. 
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where the scattering matrix components for the TM and 
TE modes are respectively given by 
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Here ..), hy{...) and h^' {. 
functions of J-th order, and k 



.) are spherical Bessel 
= ^JU^Lojc. Qs is 



the elastic part of the cross section responsible for the 
coherent scattering of light from the sample boundary. 
The absorption part Qa is responsible for diffusely scat- 
tered light via the incoherent channel and the sum Qq = 
Qs + Qa is the total cross section for the entire scat- 
tering process. Light emerging the sample via incoher- 
ent channels can be recovered in the Maxwell theory by 
considering secondary and multiply scattered waves gen- 
erated by fluctuations of This can be simulated 
by a Monte-Carlo scheme and crucially requires that the 
dipole sources generating these waves be indistinguish- 
able on a mesoscopic distance ~ A. Then secondary and 
multiply scattered waves can be simulated through the 
time dependence of the sample fluorescence when the ex- 
citation light is turned off. The analysis of the transient 
process and time dependence of the fluorescence is sim- 
plified by a diffuse approximation if the extinction length 
for the field penetration inside the sample is 3> A. 



We now turn to the quantum-posed description of the 
photon scattering problem. There photon scattering on 
an atomic system is expressed by the following relation 
between the total cross section and the T-matrix 
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where V is a quantization volume. Here we keet only the 
Rayleigh channel and assumed that the atomic system is 
described by the same ground state g before and after 
the scattering process, which includes the averaging over 
initial and sum over final Zeeman states. The initial en- 
ergy of the entire system Ei is given by Ei = Eg + fiuj 
and the incoming and outgoing photons have the same 
frequency uj' = lu. We consider below the simplest rele- 
vant example of a "two- level" atom, which has only one 
sublevel in its ground state and three Zeeman sublevels 
in its excited state, such that Fq = and F = 1 and 
in an isotropic situation the total cross section does not 
depend on the momentum direction and polarization of 
the incoming photon. 

The T-matrix is expressed by the total Hamiltonian 
H = Hq + V and by its interaction part V as 
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In the rotating wave approximation the internal resolvent 
operator contributes to (|4]) only by being projected on the 
states consisting of single atom excitation, distributed 
over the ensemble, and the vacuum state for all the field 
modes. Defining such a projector as P the projected 
resolvent 
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performs a 37V x 3iV matrix, where N is the number of 
atoms. For the dipole-type interaction between atoms 
and field this projected resolvent can be found as the 
reversed matrix of the following operator 



R-\E) ^ P [ E - Ho-VQ 



1 



QV P (7) 



E — Hq 

where the complementary projector Q = 1 — P, op- 
erating in the self-energy term, can generate only two 
types of intermediate states: a single photon + all the 
atoms in the ground state; and a single photon -I- two 
different atoms in the excited state and others are in 
the ground. For such particular projections there is the 
following important constraint on the dipole-type inter- 
action V: PVP — QVQ = 0. Due to this constraint 
the series for the reversed resolvent ([7]) is expressed by 
a finite number of terms, which would be not the case 
in a general situation The resolvent R{E) can be 

numerically calculated and, for an atomic system con- 
sisting of a macroscopic number of atoms, when A'^ 3> 1, 



3 



the microscopically calculated cross section Q can be 
compared with the Debye-Mie theory results ([2]). Both 
the microscopic and self-consistent approaches have no 
fitting parameters. Compared with the mesoscopic ap- 
proach, the essential difference is that in the microscopic 
approach any type of collective atomic excitations can 
contribute, which are not responsibly cooperative. One 
can then expect significant difference in their predictions 
for the cross section and for the fluorescence behavior. 
The time dependence of the fluorescence signal can be 
built up via Fourier expansion of the outgoing photon 
wavepacket with the S'-matrix formalism. 

Below we consider the spectral dependence of the cross 
sections, calculated in both the self-consistent mesoscopic 
model and in an exact microscopic approach. Our cen- 
tral idea is to follow how this dependence is modified 
when the dimensionless density of atoms nt)}? is varied 
from smaller to greater values. In Fig. [T]we reproduce 
such a spectral dependence when the density of atoms 
is small, noX^ = 0.02. There is excellent agreement 
between the Debye-Mie and microscopically calculated 
data and the microscopic result is insensitive to configu- 
ration averaging over random atomic distributions. This 
justifies that only cooperative modes, which allow the 
macroscopic Maxwell description, contribute to the scat- 
tering process. As a consequence the long term time 
dynamics of the fluorescence signal can be very well ap- 
proximated by a Holstein mode. Such a spectral and 
temporal behavior is typical for dilute atomic systems as 
was verified by our numerical simulations done for the 
atomic ensembles of different sizes and consisting of dif- 
ferent numbers of atoms. While varying the density to 
greater values n^)? ^ Uc ^ 0.09 the solution of the self- 
consistent equation ([T]) turns the permittivity to nega- 
tive values in a part of the spectrum where e' < and 
e" — 0. As is well known in, e.g. plasma physics ll| the 
negative permittivity can be associated with a forbid- 
den spectral zone, where the radiation cannot penetrate 
inside the system and can exist only in the form of a 
surface light-matter wave. In the limit of even higher 
density samples, the light undergoes mainly surface scat- 
tering and the absorption (incoherent/diffusion scatter- 
ing) channels suppressed. This is illustrated in Fig. [Uby 
those spectral dependencies of the cross section, which 
were calculated in the Debye-Mie model for uqA^ = 0.5. 
In contrast with similar results in Fig. [1] the contribu- 
tion of the elastic channel now dominates such that a 
smaller portion of light can penetrate the sample. As a 
consequence, it is expected that fewer atomic excitations 
should contribute to the time dependent fluorescence de- 
cay and the self-consistent model predicts an even faster 
decay than would take place in a similar dilute system. 

The situation can be differently posed if it is consid- 
ered from a microscopic point of view. Then such a dense 
atomic vapor can be relevantly described as a micro- 
cavity system built up in an environment of randomly 
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FIG. 1: The spectral dependence of the total Qo, elastic Qs, 
and absorption Qa cross sections for an atomic sample of 
size a = 25A and density noA'^ = 0.02. The configuration 
averaged microscopic result, shown as a gray solid line, is ev- 
idently the same as the data taken for particular atomic con- 
figuration, shown as a black line. The results of a mesoscopic 
self-consistent approach (red curves) reproduce the exact mi- 
croscopic spectra. Color online. 
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FIG. 2: Same as Fig[T] but for a sample size a = lOA and den- 
sity noA'^ = 0.5. The results of the self-consistent mesoscopic 
approach (red curves) perceptibly differ from the exact mi- 
croscopic spectra. In turn, the microscopic calculation given 
for a particular configuration (black curve) indicates speckle 
micro-cavity structure generated by the resolvent poles. Color 
online. 



distributed atomic scatterers. Any distribution creates 
a specific quantization problem for the incoming field, 
whose mode structure can be properly defined in terms 
of standard scattering theory. The main difficulty for the 
quantization procedure is description of the complicated 
structure of the resonance states. The resonances are de- 
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FIG. 3: Time dependent fluorescence decay for samples of 
difi'erent densities and with relatively equal optical depths. 
The dilute system has a sample radius a = 27X and density 
noA^ = 0.015, and for the dense system these characteristics 
are given by a = 5A and noA^ = 0.5. Color online. 

scribed by the resolvent poles and can be specified by 
various superpositions of atomic states, which transform 
the reversed projected resolvent ([7]) to diagonal form. 
The micro-cavity structure manifests itself by finely re- 
solved speckle dependence in the scattering spectrum and 
the scattering process becomes extremely sensitive to the 
spatial configuration of atomic scatterers, as illustrated 
in Fig. [21 Some of these resonance states have a sub- 
radiant nature and manifest themselves via significantly 
slowed long time decay of the fluorescence in compari- 
son with a classical Holstein mode. Fig. [3] illustrates the 
difference in the fluorescence decay for dilute and dense 
atomic systems initially excited by a short probe pulse. 
The external parameters for the systems are such that 
the optical depth for both the systems is nearly the same. 
For a dilute system the long term asymptote, extracted 
through exact microscopic calculations, is well described 
by a diffuse Holstein-type mode evaluated via the sample 
macroscopic characteristics. For the dense system there 
is evident deviation of the asymptotic behavior from such 
behavior. Such a deviation should be associated with the 
presence of sub-radiant resonance states in the resolvent 
poles. 

The analogy of sub-radiant states with a localization 
process, usually discussed as a multiple wave scattering 
problem, suggests evolution with density of the qualita- 
tive scattering properties of the atomic gas. The dis- 
cussed characteristics, such as the spectral dependence 
of the sample cross-section or the time-dependent fluo- 
rescence, have no direct relations with static thermody- 
namic properties of the system, which could be evaluated 
based on its partition function. However, they describe 
how the system in equilibrium responds to optical exci- 
tation near the atomic resonance transition. The atomic 
transition is dressed by cooperative dipole-field interac- 
tions and at higher densities, the disordered atomic sys- 



tem reveals a micro-cavity structure. Such a micro-cavity 
has unique properties for any particular configuration 
and its mode splitting is competitive with the mode de- 
cay rate. The modification of the dynamic response to an 
external field with density is an intrinsic property of the 
system, which may be visualized in terms of a cross over 
or a change of phase. In this sense our calculations show 
that the transformation of the system behavior from one 
of individual, independent atomic scatterers to the co- 
operatively organized micro-cavity structure mainly de- 
velops within a narrow and critical density zone. The 
configuration sensitive speckle structure of the spectral 
cross section manifests itself in those conditions when 
the self-consistent permittivity becomes negative in part 
of the spectrum. The description of the atomic system 
with the macroscopic Maxwell theory preferably yields 
surface scattering of light in this case. More precise mi- 
croscopic description shows that part of the excitation 
can penetrate the system and be converted into a long- 
time decay of the fluorescence signal. 
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